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Abstract 

We examine the effects of the additional term of the type ~ e N P N n on ^he recently proposed 
'~r [ empirical formula for the lowest excitation energy of the 2 + states in even-even nuclei. This study 
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is motivated by the fact that this term carries the favorable dependence of the valence nucleon 
numbers dictated by the N p N n scheme. We show explicitly that there is not any improvement in 
reproducing E x (2f) by including the extra N p N n term. However, our study also reveals that the 
excitation energies E x (2f), when calculated by the N p N n term alone (with the mass number A 
dependent term), are quite comparable to those calculated by the original empirical formula. 
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Recently, the role played by the valence nucleon numbers N p and N n in determining 
the lowest excitation energy of the natural parity even multipole states in even-even nuclei 
has been emphasized [l|. Over the past four decades, the valence nucleon numbers N p 
and N n have been frequently adopted in phenomenologically parameterizing various nuclear 
observables. It was Hamamoto who first pointed out that the square roots of the ratios of 
the measured and the single particle B(E2) values were proportional to the product N p N n 
J2j. Subsequently, it has been shown that a very simple pattern emerged when the nuclear 
data on the lowest excited states in the nuclei was plotted in terms of the product N p N n 
[21]. For example, the measured excitation energies E x (2f) of the lowest excited 2 + states 
in even-even nuclei make up the scattered irregular points when they are plotted against 
the mass number A (A-plot). However, they become very neatly rearranged points when 
they are plotted against the product N p Nn (NpN n -plot) [4J. A similar simplification was 
observed from the ratio E x (4f) / E x (2f) J5|, [g, |7|, the transition probability B[E2\2\ — > 



+ ) [8(, and the quadrupole deformation parameter 62 [9|. This interesting phenomenon 
has been called the N p N n scheme in the literature. In fact, the N p N n scheme has been 
extensively and successfully used for more than two decades to correlate large amounts of 
data on the collective degrees of freedom in nuclei [lOj. It has also been established that 
the N p N n scheme is manifested because the valence proton-neutron (p-n) interaction is the 
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dominant controlling factor in the development of collectivity in nuclei 
most recently the measured p-n interaction distribution in even-even nuclei over the entire 
nuclear chart has been successfully accounted for by the nuclear density functional theory 
calculations [14]. 



In a previous publication 15|, we proposed an empirical formula in the form of 



E,r. 



aA"" 1 + p 



e ~ XN p _|_ e ~ AAr « 



(1) 



for the excitation energy E x (2f) of the first 2 + states in even-even nuclei. The model 
parameters a, j3, 7, and A can be fixed easily and unambiguously by the usual least squares 
fitting technique. However, since the excitation energies span a broad range and sometimes 
the difference between the measured and calculated excitation energies can be too large, we 
take the logarithmic error Re{i), for the zth data point, of the calculated excitation energy 
E^ l (i) with respect to the measured excitation energy £^ xp (i) by [16] 



R E (i) = log 






log 



Ef(i) 



log [££*(!)]. 



(2) 



Then we minimize the dimensionless \ 2 value which is defined in terms of the logarithmic 
error by 

1 N ° I 2 

X 2 = T E\Re^\ (3) 

where N is the number of total data points considered. The empirical formula, Eq. (JTJ, has 
been found to be capable of descri bing; the essential trends of E x (2f) in even-even nuclei 



throughout the entire periodic table [15|] . We have also shown that the source, which governs 
the 2f excitation energy dependence given by Eq. ([T|) on the valence nucleon numbers, is 
the effective particle number that participates in the residual interaction from the Fermi 



level 17] . In addition, it is shown that Eq. ([I]) complies with the requirement of the N p N n 
scheme, even though it does not depend explicitly on the product N p N n [4J]. Furthermore, 
the same formula could be successfully adopted in describing the lowest excitation energy 
of the natural parity even multipole states up to 10 + in even-even nuclei [l]. 

In spite of the fact that the above empirical formula complies well with the N p N n scheme, 
there still remains some possibility that the excitation energy could depend explicitly on the 
product N p N n as dictated by the N p N n scheme, which assumes the active contribution of 
the valence p-n interaction. In this work, therefore, we want to include an additional term 
of the type ~ e~ A N P Nn (N p N n term) in Eq. (Dp) and to check how this term competes with 
other terms in reproducing the lowest 2 + excitation energy E x (2f) in even-even nuclei. The 
most general form of the proposed formula with the additional N p N n term can be written 
as: 

E x = aA~~< + (3 p e~ x ^ + ^ n e~ XnNn + p e - x ' N * Nn (4) 

where the parameters (3 and A in Eq. flTJ are split into (3 p , (3 n , and X p , A n , respectively. 
This takes into account the possibility that the contributions to the excitation energy E x 
from protons and neutrons might be different. The parameters in Eq. (j4]) are determined 
by minimizing the \ 2 value given by Eq. ([3]). Because the fitting procedure adopting the \ 2 
value in terms of the logarithmic error is quite competent, we anticipate naively that, if the 
N p N n term contributes more effectively than the separate N p and N n terms in Eq. (J3J), then 
the fitting procedure would return relatively smaller values for the parameters j3 p and (3 n 
than for (3' in Eq. (j4]). 

We perform the fitting procedure under the following five different cases, according to the 
way in which the parameters in Eq. (TJ| are treated: (I) (3 P = (3 n = (3 and X p = X n = A with 



P' = and A' = (4+0 parameters); (II) (3 P = (3 n = (3 and \ p = \ n = A (4+2 parameters); 
(III) j3' = and A' = (6+0 parameters); (IV) no restriction (6+2 parameters) and; (V) 
P P — Pn = and A p = A n = (2+2 parameters). We have distinguished the above cases by 
labeling them with two digits connected by a plus sign such as 'ra 1 + ra 2 parameters', as shown 
in the above parenthesis. The first digit denotes the number of parameters adopted from 
the first three terms in Eq. (j4"]) ? while the second digit denotes the number of parameters 
adopted from the last term of the same equation. 

The resulting parameter values for each case are listed in Tab.[U together with the corre- 
sponding dimensionless \ 2 values. Note, however, that the definition of the \ 2 value given 
by Eq. ([3]) is unrelated to one usually employed in the standard error analysis. Therefore 
the x 2 values shown in Tab. [I] should be referred for the purpose of comparison only. The 
measured excitation energies E x (2f) in this calculation are quoted from the compilation in 



Raman et al. [18j. When the N p N n term is not included in Eq. (jl]), a unique parameter set 
exists, which minimizes the \ 2 value. The cases of I (4+0 parameters) and III (6+0 param- 
eters) in Tab. [2, both correspond to this type of situation. But, we have noticed that when 
the N p N n term is present in Eq. PJ, there are two parameter sets, both of which constitute 

TABLE I: Values for the parameters in Eq. ([3]) for the excitation energy of the first 2 + state 
determined by minimizing \ 2 value under the five cases defined in the text. 

Case I II II' III IV IV V 

Parameters 4+0 4+2 x 4+2 2 6+0 6+2 x 6+2 2 2+2 



a(MeV) 


81.39 


69.65 


11.00 


68.38 


64.31 


14.45 


68.82 


7 


1.375 


1.331 


0.685 


1.342 


1.318 


0.795 


1.319 


P P (MeV) 


0.958 


0.391 


0.922 


0.826 


0.399 


0.796 




/3n(MeV) 


0.958 


0.391 


0.922 


1.167 


0.754 


1.141 




Xp 


0.344 


oc 


0.310 


0.419 


0.750 


0.359 




*n 


0.344 


oc 


0.310 


0.284 


0.374 


0.284 




/?'(MeV) 




0.795 


-0.311 




0.510 


-0.216 


0.964 


A' x 10 4 




233 


9 




215 


11 


271 



2 0.151 0.128 0.126 0.126 0.117 0.109 0.140 



X 



two local minima for the x 2 value. These two parameter sets are distinguished in Tab. [J by 
adding a subscript to the second digit (e.g., 4+2i and 4+2 2 ). The cases of II (4+2i and 
4+2 2 parameters) and IV (6+2! and 6+2 2 parameters) correspond to this type of situation. 
However, when only the N p N n term is present in Eq. (|5j), such as in case V in Tab.[U there 
is a recurrence of a unique parameter set. 

Fig. [1] presents the results for the first 2 + excitation energies E x (2f) in even-even nuclei 
obtained by the four-parameter empirical formula, Eq. ([T|). The results obtained by the 
six-parameter empirical formula, Eq. (j4"|), are presented in Fig.[2j In both figures, the points 
connected by solid lines represent nuclei that belong to the same atomic number (the isotopic 
chains). The dashed lines of the two figures are drawn by taking only the first term aA~"< 
of the empirical formula. The calculated excitation energies are to be compared with the 
measured data shown in Fig.[T](a). 

In Figs.[T](b), 1(c), and 1(d), the A-plots of the excitation energies E x {2±) are shown, 
which are calculated by adopting the parameter set of the 4+0 parameter for Fig. 1(b), 
the 4+2i parameter for Fig. 1(c), and the 4+2 2 parameter for Fig. 1(d) as given in Tab.[H 
As has been discussed in the previous publications [l|, ll5(, we find from Fig. 1(b) that the 
four-parameter empirical formula, Eq. (JT]), describes reasonably well the essential trends of 
the lowest 2 + excitation energies in even-even nuclei throughout the entire periodic table. 
We also find, by comparing Figs. 1(c) and 1(d) with 1(b), that the N p N n term added to 
the four-parameter empirical formula, does not have a significant influence on the overall 
shape of the graphs. However, we also notice that there are a few small variations brought 
about by the N p N n term. First of all, we find that the N p N n term can contribute either 
additively (as in Fig. 1(c)) or subtractively (as in Fig. 1(d)) to the original formula. This 
fact can be reconfirmed by the opposite sign of the parameter /3' for the cases II and II' in 
Tab.[Tl It is of notable interest that the dashed curve in Fig. 1(c) (where the N p N n term adds 
to the remaining terms), follows the bottom contour line of the excitation energies, while 
the dashed curve in Fig. 1(d) (where the N p N n term subtracts from the remaining terms), 
passes through the middle of the excitation energies. Even though the x 2 value for case II' 
(4+2 2 parameters) is lower than it is for case II (4+2i parameters) by only 0.002, we find 
that the results obtained by case II reproduce the measured excitation energies better than 
the results obtained by case II' in the mid-shell region. However, it is interesting to find 
that, in case II, there is no contribution from the second and third terms in Eq. (j4"l), except 
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FIG. 1: Excitation energies of the first 2 + states in even-even nuclei. The points are connected 



by solid lines along the isotopic chains. Part (a) shows the measured excitation energies Ha] while 
parts (b), (c), and (d) show those calculated by our empirical formula given by Eq. Q with the 
parameter set given by 4+0, 4+2i, and 4+22 parameters in Tab.[U respectively. The dashed lines 
in parts (b), (c), and (d) are drawn by taking only the first term a A'" 1 of Eq. ([4]). 
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FIG. 2: Same as Fig. [T] but for parts (a), (b), (c), and (d) show the excitation energies calculated 
by our empirical formula, Eq. §%§ with the parameter set given by 6+0, 6+2i, 6+22, and 2+2 
parameters in Tab.[U respectively. 



for N p = or N n = 0, as a result of the fitted parameter values X p = X n = oo. Consequently, 
an almost straight line represents the A-plots of the excitation energies E x (2f) for N p = 0, 
while the excitation energies E x (2f) for N p = N n = nuclei are represented by dots above 
the previous straight line in Fig. 1(c), which does not concur with the shape given by the 
data. This eccentric behavior occurs because, under the condition of j3 p = f3 n and \ p = A at, 
the last N p N n term can play the role of both the second and third terms in Eq. (j4j). 

By comparing case I with case II in Tab. [J we learn that, when the N p N n term is added 
to the four-parameter empirical formula, the \ 2 value is lowered by approximately 15%, 
from 0.151 to 0.128. Almost the same amount of reduction in the \ 2 value can be obtained 
if the four-parameter formula is extended to the six-parameter one. This can be seen from 
the x 2 values of case I and case III in Tab.[H However, by adding the N p N n term to the six- 
parameter empirical formula (compare case III with case IV), a further reduction in \ 2 value 
can be achieved, which amounts to only approximately 6%. In Figs. 2(a), 2(b), and 2(c), 
we show the A-plots of the excitation energies E x (2f). These are calculated by adopting 
the parameter set of the 6+0 parameters for Fig. 2(a), the 6+2i parameters for Fig. 2(b), 
and the 6+22 parameters for Fig. 2(c), as given in Tab.Hl The three graphs calculated by 
the six-parameter empirical formula in Fig. [2] are basically same as the corresponding three 
graphs calculated by the four-parameter counterparts in Fig. [1] except for the graph shown 
in Fig. 2(b). For case IV, where the four parameters j3 p , j3 n , X p , and A n are separately fitted, 
there are suitable contributions from the second and third terms in Eq. (Ill) and the eccentric 
behavior shown in Fig. 1(c) disappears, as shown in Fig. 2(b). 

By reviewing Figs. Q] and El it becomes apparent that the role played by the supple- 
mentary N p N n term is insignificant when it is added to the four or six parameter empirical 
formula. However, we present a very interesting results in Fig.El(d) where we show the 
A-plot of the excitation energies E x (2f) calculated by only the N p N n term together with 
the first a A -1 term in Eq. (j4j). In this figure, we adopt the 2+2 parameter set which is 
listed as case V in Tab.fl] It is quite surprising to find that the graph obtained by the 2+2 
parameter set is highly comparable to that obtained by the 6+2i parameter set, with the 
exception of the excitation energies in nuclei that belong to the closed proton and/or the 
neutron major shell. Furthermore, since the \ 2 value should decrease monotonically as the 
number of parameters increases, the results of the 2+2 parameter set, where x 2 = 0.140, 
could be regarded better than, or at least as good as, the results of the 4+0 parameter set, 

8 



where \ 2 = 0.151. 

In order to provide a closer confirmation of the relative performance of the various cases 
presented in Figs.[TJ and [2J we choose a number of individual isotopic chains and remake 
the A-plots of them in Figs. [3] and HI In Fig.[31 we display arbitrarily chosen four isotopic 
chains, Z = 30, Z = 36, Z = 42, and Z = 48, which include nuclei that belong to the 
proton major shell between Z = 30 and 50. The left panels of Fig.[3] show the measured 
excitation energies E x (2+) (solid squares) as well as those calculated by adopting the 4+0 
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FIG. 3: Excitation energies of the first 2 + states in nuclei which belong to the isotopic chains 
Z = 30, Z = 36, Z = 42, and Z = 48. The measured excitation energies (solid squares) as 
well as the calculated ones by adopting the 4+0 (solid triangles), 4+2i (solid circles), 2+2 (empty 
squares), 6+0 (empty triangles), and 6+2i (empty circles) parameter sets are shown. 



(solid triangles) and 4+2i (solid circles) parameter sets. In addition, the right panels of the 
same figure show the excitation energies E x (2f), calculated by adopting the 2+2 (empty 
squares), 6+0 (empty triangles), and 6+2x (empty circles) parameter sets. These graphs 
confirm that the addition of the N p N n term has a minimal effect in both empirical formulae, 
Eq. (pQ) and Eq. (Pi]) . Furthermore, it is interesting to discover that the excitation energies 
E x (2f) calculated by only the N p N n term with the mass number dependent aA~^ term (i.e., 
by the 2+2 parameter set), are as satisfactory as those calculated by any other parameter 
sets. Similarly in Fig.HJ we display again arbitrarily chosen four isotopic chains, such as 
Z = 62, Z = 72, Z = 78, and Z = 82, which include nuclei that belong to the proton major 
shell between Z = 52 and 82. This figure also shows that the excitation energies E x (2i), 
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calculated by the parameter sets with the additional N p N n term, are found to be more or 
less the same as those that were calculated by the parameter sets without the N p N n term. 
While the differences between these graphs are small, the case with the best results is the 
6+0 parameter set. Fig. H] also shows that the results obtained by the 2+2 parameter set are 
very similar to those obtained by any other parameter sets, with the exception of Z = 82, 
which is the proton closed shell. 

Finally, we plot the excitation energies E x (2f) shown in Figs.H] and [2] again m Fig.[5]but 
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FIG. 5: The -/VpiV n -plot versions of the ^4-plots shown in Fig.Q]and[2j The excitation energies are 
expressed by the following different symbols according to which major shells they belong to: solid 
squares (Z = 2 ~ 28), empty circles (Z = 30 ~ 50), solid circles (Z = 52 ~ 82), and empty squares 
(Z = 84 ~ 100). 
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this time we make the iV p iV n -plots instead of the A-plots. Of course, the graphs in Fig. [5] 
are drawn with exactly the same set of plotted points as those that were used in Figs. CD 
and [21 The three left panels of Fig.Oshow the iVpJV n -plot versions of the A-plots that were 
shown in Fig.[T](a) (data),[TJb) (4+0), and[TJc) (4+2!), while the three right panels show the 
JVpiV n -plot versions of the A-plots shown in Figs.^d) (2+2),[2](a) (6+0), and[^b) (6+2i). 
In these iV^A^- plots, the excitation energies E x (2f) are expressed by the following different 
symbols, according to which major shell they belong to: solid squares (Z = 2 ~ 28); empty 
circles (Z = 30 ~ 50); solid circles (Z = 52 ~ 82); and empty squares (Z = 84 ~ 100). We 
can observe from these N p N n -plots (as with the A-plots), that the overall characteristics of 
the data are reasonably reproduced by our empirical formula. We also find from the N p N n 
plots that the N p N n term added to the empirical formula does not alter the essential shape 
of the graph. However, we can still recognize the effect brought by the addition of the N p N n 
term to the empirical formula by comparing the lowest two left panels or the lowest two 
right panels. The width of the plotted points in the N p N n -p\ot becomes narrow and the 
plotted points which belong to the different major shells become further separated when 
the N p N n term is present, than when the N p N n term is absent. However, by comparing 
calculated results with the data shown in the top left panel of Fig. [5], it is evident that both 
of these effects, caused by the additional N p N n term, deteriorate the fit to the data. Just as 
pointed out regarding Fig. [2], the N p N n plot drawn by the 2+2 parameter set in Fig.[5]is very 
similar to the iVpiN^-plot drawn by any other parameter sets. Furthermore, it is identical to 
the N p N n -p\ot drawn by the 4+2i parameter set. This is because, as discussed previously, 
the second and third terms of Eq. (T4J) do not contribute in any way to E x (2i), when this 
parameter set is adopted. 

In summary, we wanted to double check any possibility that the first 2 + excitation en- 
ergy could depend explicitly on the product N p N n as expected by the N p N n scheme. We, 
therefore, included the N p N n term of the type ~ e~ A N P Nn as in Eq. (T4J) and investigated how 
this term competed with other terms in reproducing the first 2 + excitation energy E x (2f) 
in even-even nuclei. In fact, we anticipated naively that if the N p N n term contributed more 
effectively than the separate N p and N n terms in Eq. (T4]), then our fitting procedure would 
return much smaller values for the parameters (5 P and (5 n than f3' in Eq. (T4]). But contrary 
to our anticipation, our results show that the role played by the N p N n term is insignificant 
when it is added to the former empirical formula as evidenced by the A-plots displayed in 
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Figs. [I] and [2j However, we also end up with an unexpected but very interesting result. The 
A-plot drawn by the 2+2 parameter set (calculated with only the N p N n term and the ctA~ 7 
term) is quite similar to those drawn by any other parameter set. But the result calculated 
by the 2+2 parameter set has the following two observational flaws. One is that it could not 
reproduce the excitation energies E x (2f) in nuclei that belong to the closed proton and/or 
neutron major shell where N p = and/or N n = 0. The other is that the N p N n -plot of the 
2+2 parameter set is too narrow in order to match the N p N n -plot of the data. However, one 
must not forget the fact that the x 2 value for the 2+2 parameter set is quite comparable to 
that for the 4+0 parameter set. In fact, it is fair to say that the present results alone could 
not tell which one, the N p cross N n term or the N p plus N n term, would finally represent 
the lowest excitation energy systematics in nuclei. A further study to clear up this situation 
is called for. 
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